1 Introduction 

The ability to prepare arrays with made to measure electronic properties 
has generated much recent activity 0]. This was made possible by the de- 
velopment of synthetic methods for the preparation of, so called, quantum 
dots. For our purpose these dots act like designer atoms. The properties of 
individual dots can be designed by the choice of composition and size. Size 
is a particular variable we wish to emphasize. Quantum dots are aggregates 
of tens or hundreds of atoms (or molecules) and so are nanoscale building 
blocks. The size of the dots is however small enough so that, since the high- 
est lying electron(s) are confined to the volume of the dot, the electronic 
energy levels are still discrete. For a chemist, it is convenient to think of the 
higher-most electrons of each dot as the valence electron of an atom|I. 

These dots can self assemble into a planar array. For, e.g., Ag nanodots, 
the packing is hexagonal. The reason why adjacent dots in such an array may 
have an unequal charge distribution is related to their nano-scale size. Unlike 
the smaller ordinary atoms, a quantum dot does not require a high energy 
cost to induce an electron transfer from a neighboring dot. This is techni- 
cally stated by saying that the charging energy of a dot is atypically low. 
Electrostatic considerations suggest that the charging energy scales inversely 
with size and this has been experimentally born out. Metallic quantum dots 
have charging energy of well below an eV while atoms have charging energies 
of at least several eV's. 

The current attention on architecture with quantum dots building blocks 
is partly motivated by potential applications of such nanoscale quantum de- 
vices. In this paper we explore how corresponding computations can be 
carried out. The technical problem is that the Coulombic repulsion between 
two electrons (of opposite spins) that occupy the same dot cannot be de- 
scribed in a one electron approximation. It requires allowing for correlation 
of electrons. The simplest model Hamiltonian that can describe the Coulom- 
bic (or charging energy) effects is due to Hubbard 0. This paper is a study 
of exploring a computationally tractable method for solving for the states of 
the Hubbard Hamiltonian for a lattice. The development makes essential use 
of the nature of the intended applications. At low temperatures, only those 
states of the array that are quite close to the ground state are important ||. 

The Hubbard model is the simplest model one can use to study the Mott 
metal-insulator transition, i.e. the transition from a metallic to an insulat- 
ing phase driven by correlations between electrons [|5[]. Despite the extensive 
efforts in studying the Hubbard model || ||, [7| there are still difficulties in 
extending the one-dimensional results to higher- dimensions. The difficulties 
mainly come from the electron correlations inherent in this model, which 
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becomes more complicated in two and three-dimensions. Several methods 
have been developed to treat these problems such as exact diagonalization 
met hods H|, quantum Monte-Carlo simulations and a variety of approx- 
imation techniques like mean-field theory |TD|, Green's function decoupling 
schemes []Ll|], functional integral formulations fl2|], variational approach [T3f 



and perturbation expansion ||14j| . It is already known that the numerical 
methods suffer from the finite-size effect and perturbation calculations are 
not useful when dealing with the intermediate-to-strong correlations. In this 
context, using the real-space renormalization group method seem to be a 
promising direction because of its non-perturbative nature. In fact, this 
method has been applied to quantum lattice systems flB]] . For example, the 
density matrix renormalization group \W\ has provided a revolutionary way 
to obtain reliable results for the one-dimensional Hubbard model and the 
numerical renormalization group for interacting finite Fermi systems [|T7|] and 
for excitations in atoms|18|]. To apply the renormalization group method 
to higher dimensions, the simplest way is to use the block renormalization 
group method(BRG)||19|. The BRG method has already been applied to in- 



teracting electrons in two-dimensional systems for both bipartite lattices [20 



and more recently for nonpartite lattices such as the triangular one pi 



In this paper, we present the real-space block renormalization group 
equations for fermion systems on a triangular lattice with hexagonal blocks. 
By using operator transformations, we obtain the conditions that make the 
renormalized Hamiltonian in the bare Hubbard model have the same struc- 
ture as the old one, i.e. the renormalization group equations in the pa- 
rameters space with constant dimension will remain the same, without any 
proliferation of couplings. Finally, using numerical calculations we show that 
the metal-insulator transition occurs at the critical value of U/t ~ 12.5. 



2 Real-space renormalization group 

In this section, we present a general computational procedure for applying 
the BRG method for interacting fermions on a lattice using the Hubbard 
Hamiltonian. As an example we will take the hexagonal block on a triangular 
lattice. The Hubbard Hamiltonian can be written as follows, 



H = -t J2 [4c c jT + H.c.\ + uj2 n n n n -/*D n <t+ n *i)' W 

<i,j>,a i i 

where t is the nearest-neighbor hopping term, U is the local repulsive 
interaction and /i is the chemical potential. cf a (c ia ) creates(annihilates) an 
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electron with spin a in a Wannier orbital located at site i; the corresponding 
number operator is n ia = cf a c io - and < > denotes the nearest-neighbor pairs. 
H.c. denotes the Hermitian conjugate. 

The essence of the BRG method is to map the above many-particle Hamil- 
tonian on a lattice to a new one with fewer degrees of freedom and with the 
same low-lying energy levels [2^|. Then the mapping is repeated leading to 



a final Hamiltonian for which an exact solution can be obtained. The pro- 
cedure can be divided into three steps: First divide the N— site lattice into 
appropriate n s — site blocks labeled by p(p = 1, 2, N/n s ) and separate the 
Hamiltonian H into intrablock part Hb and interblock Hjb 

H = H b + H ib = J2H p + E V p j, (2) 
v <p,p'> 

where 



H P = -t E [ c %) a c jV>o + H - C \ + U^2n i ( P )^n i ( P ) l - A^E^^t + 

<i(p) j(p)> i(p) i(p) 

(3) 

and 

V P y = ~t E [ct^^'v* + H.C.] (4) 
<j(p,i>)j(p',i))> 

in which denotes the ith site on the pth block and i^'^ denotes the 
border site of the block p. 

The second step is to solve H p exactly to get the eigenvalues E p . and eigen- 
functions = 1, 2, 4™ s ). Then we can build the eigenfunctions of Hb by 
direct multiplication of which can be written as |\1/b(^i^2, ■■■,iN/n s ) >— 

> \<5>2i 2 > -\®N/njN/n s > {h,h,- & {1, 2, . .. , 4" fl }) . 

The last step is to treat each block as one site on a new lattice and the 
correlations between blocks as hopping interactions. The original Hilbert 
space has four states per site. By following the above procedure one obtains 
an equivalent Hamiltonian with (zps^M _ 4^ degrees of freedom, which is 
the same as the original Hamiltonian. But in the realistic case, if we only care 
about the properties related to some special energy levels of the system it is 
not necessary to keep all the states for a block to obtain the new Hamiltonian. 



For example, when studying the metal- insulator-transition [23|, we may only 
need to consider the ground state and the first excited state energies. 

The above scheme is a general procedure for applying the BRG method. 
In order to make the new Hamiltonian more tractable, it is desirable to make 
it have the same structure as the original one, i.e. the reduction in size should 
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not be accompanied by a proliferation of new couplings. Then we can use the 
iteration procedures to solve the model. To achieve this goal it is necessary 
to keep only 4 states in step 2, which can be understood from the following 
renormalized intrasite Hamiltonian. The four selected states are taken to be 



Qpl >= 


|o>;, 




(5) 


%2 > = 


J+ In > / _ 1 


n>;, 


(6) 


%3 > = 


St i u 


r>;, 


(7) 


$ p4 > = 


c 'p + i\°>' p = 




(8) 



where c' pa (c' pa ) is the creation(annihilation)operator of the block state 
\a >' p and their corresponding energies are Ei (i = 1, 2, 3, 4). 

Our next task is to rewrite the old Hamiltonian H = Hb + Hjb in the 
space spanned by the truncated basis 



H' 



E 



I -^Truncated ^ ^ -^Truncated 



„, -Truncated -Truncated 

H\ \$/ B ><^b 



rr~\ , , —Truncated 

^fl runcated-qj ^ 



where the truncated basis is given by 



(9) 



I* 



Truncated i 



(i 1: i 2 , ...,lN/n.) >= |$l*i > \®2i 2 > —\&N/n.iN/n. > (»1,«2, --■ € {1,2,3,4}) 

(10) 

In order to avoid proliferation of additional couplings in H', the four states 
kept from the block cannot be arbitrarily chosen. Some definite conditions 
must be satisfied in order to make H' have the same structure as H. Substi- 
tuting H into H' and using the product of different operators (see Table I) 
we can get the expression for H p 



h p = |o>;s 1 <oi;+in>;^2<iTi;+u>;^4<ii;+iT>;^3<Ti; 

= E 1 + (E 3 - E 1 )n' pA + (E 4 - E 1 )n' pd + [E 1 + E 2 - E 3 - E^n'Hl) 

Note that by keeping only four states from the block states in the beginning 
gives no other extra couplings in the new Hamiltonian. 

Table I: The product of different operator transformations* 





< or 


<t r 


<i r 


<n r 


|0>' 


1 — n'l — n\ + n\n\ 


c T - n^c T 


c i - n 'A 




IT>' 


c T - c T n[ 






n T c l 


II>' 


'+ _ r '+ / 

c l c l n T 


j+ c i 


n\ - n\n\ 


n 'l c \ 




'+ >+ 


~ n 1 c l 
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*In this table the product reads |0 >'< 0|' = 1 — n' T — n'^ + n'+n'^, etc. 

Comparing the above intrasite Hamiltonian with Eq.(l), we get the next 
conditions in order to copy the intrasite structure of the old Hamiltonian, i.e. 
E 3 = £4. Because of the additional vacuum energy E\ in the new Hamiltonian 
we rewrite the intrasite part of Eq.(l) as 

H B = UJ2 n i\ n H - vJ2( n n + n n) + K J2 I i ( 12 ) 

i i i 

where we introduce another parameter K to the original system and is 
a unit operator. The new intrasite Hamiltonian is given by 

H' B = (E 1 + E 2 - 2E 3 ) £ n' p1 n' pl - (E, - E 3 ) ^(n^ + n'J + E, ^ /„. (13) 

v p p 

Then the renormalized parameters U, \x and K can be obtained from the 
following relations 



U' = E l + E 2 - 2E 3 (14) 
li' = E 1 -E 3 (15) 
K' = Ex (16) 

in which Ei, E 2 and E 3 are functions of the old parameters t, U, fi, K. 
For the half-filled case, fi — U/2, H B can be expressed as 

He = UJ2(\- ^){\ ~ mi) +Kj2h (17) 

i i 

with the initial value of K = — j. By using the particle- hole symmetry, 
Ei = E 2 , the renormalization group equations for U and K take the form 



U' = 2(E 1 -E 3 ) (18) 
K' = (Ei + E 3 )/2 (19) 

To illustrate this procedure, let us consider the triangular lattice with 
hexagonal blocks as shown in figure 1. For this non-bipartite lattice the 
interaction between blocks can be written as 
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Figure 1: Schematic diagram of the triangular lattice with hexagonal blocks. 
Only two neighboring blocks p and p' are drawn here. The dotted lines 
represent the interblock interactions and solid line intrablock ones. 



Kp' = (-*) E >< QpiA^wJ^ >< "VJ] x W%'i2 >< %'i 2 \c 5(p/)a \% % >< % % 

+ [\%h > < %h\ct(v)a\®Pi'i >< VJ1 X tl^ia >< %'i 2 \c 4 ( P ') a \% % >< %%\] 

+ [\% h > < $ pil |c+ (p)(T |$ K >< VJ1 X [l $ P'i2 >< %'i 2 \c 5 ( P ') a \% % >< ®P%\] + #C} (20 

Since we would like to keep V pp > of the form 

V PP ' = {-t')Y.K^ + H.Cl (21) 

<J 

we use the product transformation in Table I to simplify Eq. (20), 

V PP> = 2 {< < C iWl° >p +[< - ff > ^cjrfj - CT >P - < ^IpC^J > P K— 
a,<i,j> 

X{ < 0|p,C j(p O CT |(T >p, +[< -(Tl^C^ak, -° IT>p' - < ^Ip'C^'U >p'K'-J C p'a + 

( < ij >=< 1,5 >,< 2,4 >,< 2,5 >). (22) 

It can be easily seen now that in order to make all the extra couplings vanish, 
it is necessary to make further restrictions upon the selected states, 
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< - (J , (J \p c t(p)a\p- a >'p =< a \'p c t)J° >p> ( 23 ) 

< -vljCjWaW, -° >' P '=< %Cj{v')M >' p , (24) 

Using similar calculations to the other neighboring interactions of the block, 
we can finally obtain the following conditions, 

< - ff > vl'pcfwjp - ° >'p=< vlpcfwjQ >' P = x ( 25 ) 
for all the border sites on the block. Then the new hopping term becomes, 

V pp/ = u\^c' p %^ (26) 

where v represents the number of couplings between neighboring blocks. In 
figure 1 v — 3. The last renormalization group equation is readily obtained 

t' = v\H. (27) 

Up to now, we have given a general discussion of the conditions under 
which no proliferation of couplings results from the application of the BRG 
method to non-partite lattice. Because on the border of non-partite lattice 
block, there is only one type of sites, the above procedures can be extended 
to other lattices with different dimensions or blocks without much difficulty 



3 State Selections 

After deriving the conditions for the renormalization group equations, the 
next task is to select states that satisfy these conditions. At this stage the 
symmetry properties of the lattice play an important role. From Eqs.(23) 
and (24), it can be easily seen that if we assume the particle number in the 
state |0 >' to be N e -1, then in |f>', | |>' and | f j>', there should be N e , N e , 
N e + 1 particles respectively. Moreover if the spin in |0 >' is S z , the spins for 
|t>',| |>'and| tl>', should be S z + l/2, S z -l/2 and S z . The total electron 
number N e and the spin S z for each block are good quantum numbers since 
their corresponding operators commute with the Hubbard Hamiltonian. So 
when we diagonalize the Hubbard Hamiltonian of the selected block, we keep 
iV e and S z fixed to be (iV e -l, S z ), (N e , S g +l/2), (iV e , S z -l/2) and (N e +1, 
S z ) respectively. Thus we obtain four groups of eigenenergies and eigenstates 
corresponding to the above quantum numbers. From each group, we select 
the lowest-energy state to form the final required four states. It should 
be mentioned that the lowest-energy state has to be selected according to 



7 



definite special symmetry considerations, which shall be discussed in the 
next paragraph. In order to obtain the insulating to conducting gap, which 
is defined to be the energy difference between extracting one electron from 
the system and adding one electron to it, N e is selected to be equal to N s . 
For S z , we choose it to be zero so as to make the block have the same spin 
property as the one-site. So now the renormalized lattice will be composed 
of N/n s renormalized "sites" with N/n s "particles". 

Instead of forcing the above conserved quantities upon the selected states 
in analogy to the one site properties in a consistent way, here we get them 
directly from the no-coupling-proliferation conditions. A does not depend on 
a in Eq. (25), this can be guaranteed by the particle- hole symmetry, which 
means that only in half-filled lattice can the renormalized Hamiltonian have 
exactly the same form as the original one ||24|| . Moreover, the irrelevance 
of A with respect to the border site can be shown by requiring the 
selected states to belong to the same irreducible representation of the spatial 
group of the lattice. For the triangular lattice with hexagonal blocks, the 



Hamiltonian is invariant under C% v 25 . So if we choose the same one- 



dimensional irreducible representation of the group Cq v for |0 >' , ||>', | !>' 
and | T!>'j the conditions can be satisfied. 



4 Numerical Calculations and discussions 

The insulating to conducting gap A g [[F] and the ground state energy per site 
Eg for the half-filled Hubbard model on a triangular lattice with hexagonal 
blocks can be written as 

A g = lim n ^ 00 U {n) (28) 

and 

K (n) + t/M jj 

E g = lim r ^ 0B tT^ + 2 ( 29 ) 

where n denotes the number of iterations in the renormalization equations. 
As n — > oo, the lattice containing 7 n sites will be renormalized to one 'site', 
which can be solved exactly to obatin the above equations. U/2 is added in 
Eq.(29) in order to compensate the chemical potential energy substracted in 
the original Hamiltonian. The metal-insulator-transition for this lattice can 
be examined by considering the energy gap A as a function of (U/t). If the 
insulating gap disappears the system exhibits metallic behavior. Otherwise 
it is insulating. 



8 



O 4 8 12 16 20 

u/t 



Figure 2: C//t dependence of (a) the insulating gap and (b) the ground-state 
energy per site. 



Figure 2 shows the numerical results for the scaled gap A 9 /t as a function 
of (U/t). There is a first-order phase transition at [U/t) c pa 12.5. This 
finding is quite different from the case of half-filled square lattice where an 
insulating gap exists for arbitrarily small values of U/t. This is because there 
is a perfect nesting of the Fermi surface on a square lattice, which makes the 
model unstable toward antiferromagnetism as soon as a nonzero U is turned 
on, driving the system to an insulating state. 

On the triangular lattice, for the lack of perfect nesting, metal-insulator- 
transition takes place at a finite value of U/t. But because there is no 
exact solution to the two-dimensional Hubbard model, there are still con- 
troversies upon the exact value for this critical point. The situation be- 
comes more complicated and subtle because of the frustrations inherent in 
the triangular lattice, which may induce a nontrivial competitions among 
different magnetic phases. In the Hartree-Fock calculations by Krishna- 
murthy et al.|EB[ |27| they found that for small U/t the system is a paramag- 
netic metal which turns to a metal with incommensurate spiral spin-density 
wave at U c \/t = 3.97. Two successive first-order phase transitions occur at 
U C 2/t = 4.45, a semi-metallic linear spin-density wave is stabilized and a first- 
order metal-insulator-transition to an anti-ferromagnetic insulator occurs at 
Uci/t = 5.27. Capone et al.|22| obtained qualitatively similar phase transi- 
tions by using the Kotliar-Ruckenstein slave-boson technique. That is, the 
weak-coupling paramagnetic metal continuously evolves into a spiral metal at 
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U c i/t = 6.68, which crosses the linearly polarized spin-density-wave ground 
state at V c %jt = 6.84. The latter phase undergo a further first-order transi- 
tion toward an anti-ferromagnetic insulator at U^/t = 7.68. The exact diag- 
onalization results exhibit a first-order transition between the paramagnetic 
metal and the anti-ferromagnetic insulator at U/t = 12.07, without interme- 



diate "exotic" phases [29] . Our results for the critical value of (U/t) ~ 12.5 
for the metal-insulator-transition and the ground state energy as a function 
of U/t are in agreement with results obtained by the exact diagonalization 
method ^] . 

In summary, we find that there is a first-order phase metal-insulator- 
transition, at (U/t) c = 12.5 in full agreement with the exact diagonalization 
result U/t = 12.07. 



5 Concluding Remarks 

A practical renormalization method for an array described by a Hubbard 
Hamiltonian has been presented. Specifically, we could find conditions under 
which the renormalized Hubbard Hamiltonian on a nonpartite lattice has 
no proliferation of coupling. The method is especially suited for getting the 
lowest excited states. An application demonstrating the identification of a 
first order metal-insulator transition has been discussed in detail. In future 
work we intend to build upon the present results in two directions. First, 
to apply them to the computation of the transport properties of lattices of 
quantum dots. We also hope to be able to incorporate into the method 
one essential aspect of quantum dots that, so far, is not included. This 
is the (small but finite) variability of the chemical potential. This results 
from the size fluctuations of dots as they are prepared in the laboratory. 
In real-space BRG method, how to truncate the states for one block is of 
central importance to the accuracy of the final results. In this paper, four 
states for each block are selected from the requirement of non-proliferations 
of parameters. This greatly simply the theoretical handling of this problem 
while keeping the main related physics. But for two-dimensional system, how 
to judge the convergence when more states are kept is still an open problem 
in this direction. More work is under way. 
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Figure 3: U/t dependence of (a) the insulating gap and (b) the ground-state 
energy per site. 
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